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ABSTRACT 

Massive stars influence their surroundings through radiation, winds, and supernova explosions far 
out of proportion to their small numbers. However, the physical processes that initiate and govern the 
birth of massive stars remain poorly understood. Two widely discussed models are monolithic collapse 
of molecular cloud cores and competitive accretion. To learn more about massive star formation, we 
perform and analyze simulations of the collapse of rotating, massive, cloud cores including radiative 
heating by both non-ionizing and ionizing radiation using the FLASH adaptive mesh refinement code. 
These simulations show fragmentation from gravitational instability in the enormously dense accretion 
flows required to build up massive stars. Secondary stars form rapidly in these flows and accrete mass 
that would have otherwise been consumed by the massive star in the center, in a process that we term 
fragmentation-induced starvation. This explains why massive stars are usually found as members of 
high-order stellar systems that themselves belong to large clusters containing stars of all masses. The 
radiative heating does not prevent fragmentation, but does lead to a higher Jeans mass, resulting in 
fewer and more massive stars than would form without the heating. This mechanism reproduces the 
observed relation between the total stellar mass in the cluster and the mass of the largest star. It 
predicts strong clumping and filamentary structure in the center of collapsing cores, as has recently 
been observed. We speculate that a similar mechanism will act during primordial star formation. 



1. INTRODUCTION 

Understanding the formation of massive stars is of piv- 
otal importance in modern astrophysics. Massive stars 
are rare and short lived. However, they are also very 
bright and allow us to reach out to the far ends of the 
universe. For example, th e most distant galaxie s in the 
Hubble Ultra Deep Field ( jBeckwith et al.||2006| are all 
characterized by vigorous high-mass star formation. Un- 
derstanding the origin of massive stars, at present and at 
early times, therefore is a prerequisite to understanding 
cosmic history. In our own Milky Way, high-mass stars 
contribute the bulk of the UV radiation field. This ra- 
diation can ionize hydrogen and dissolve molecules. Ex- 
panding H II regions, bubbles of ionized gas surround- 
ing massive stars, have been identified as an important 
source of interstellar turbulence. The complex interplay 
between gas dynamics and radiation thus is a key ele- 
ment of the evolution of the interstellar medium (ISM) 
and the Galactic matter cycle. Furthermore, stars are 
the primary source of chemical elements heavier than the 
hydrogen, helium, and lithium that were produced in the 
Big Bang. Again, massive stars have a disproportionally 
large share in the enrichment history of the Galaxy and 
the universe as a whole. It is the violent supernova explo- 
sions associated with the death of high-mass stars that 
contribute most of the metals. Because they insert large 
amounts of energy and momentum input, supernovae at 
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the same time strongly stir up the ISM and provide for 
the effective mixing of the newly bread elements. 

Despite their importance, the physical processes that 
initiate and control the build-up of massive stars are 
still not wel l understood a.nd subject of intense debate 
Mac Low & Klcssen 2004} [Zinnecker fc Yorke|2007| [Mc^ 



i 



ee fc: Ostriker||2007i ) . ^Because their tormation time is 
short, of order oFTO^yr, and because they form deeply 
embedded in massive cloud cores, very little is known 
about the initial and environmental conditions of high- 
mass stellar birth. In general high-mass star forming re- 
gions are characterized by more extreme physical condi- 
tions than their low-mass counterparts, containing cores 
of size, mass, and velocity dispersion roughly an order 
of magnit ude larger tha n those of cores in low-mass re- 
Jijinaet al.|19 99; Gar ay fc: Lizano|1999 Kurtz 
' 'Beuther et al.|2007j ^Motte et al.f2008[ ) Typ- 



gions (e.g. 
et al.|2000 



ical sizes of cluster-tormmg clumps are ~ 1 pc, they have 
mean densities of n '-^ 10^ cm~'^, masses of ~ 10^ Mq and 
above, and velocity dispersions ranging between 1.5 and 
4kms~^. Whenever observed with high resolution, these 
clumps break up in even denser cores, that are believed 
to be the immediate precursors of single or gravitation- 
ally bound multiple massive protostars. 

Massive stars usually form as mem bers of multiple stel- 
lar sys tems (H o fc Haschick||1981[ |Lada||2006 Zinnecker 
fc: Yorke 2007J which themselves are parts of larger clus 



ters JLada fc Lada||2003[ |de Wit et ~S'. 
1997 ). This tact adds additional cha 



2004 Testi et al. 



lenges to the in- 
terpretation of observational data from high-mass star 
forming regions as it is difficult to disentangle mutual 
dynamical interactions from the influence of individual 
stars (e.g. Goto et al. 2006 Linz et al. 20051. Further- 



more, high-mass stars reach the main sequence while still 
accreting. Their Kelvin-Helmholtz pre-main sequence 
contraction time is considerably shorter than their ac- 



cretion time. Once a star has reached a mass of about 
10 M© its spectrum becomes UV dominated and it be- 
gins to ionize its environment. This means that accretion 
as well as ionizing and n on-ionizing radiation needs to 
be con sidered in concert ( |Keto||2002[ |2003[ |2007| [Peters 
et al. |2010ab). It has been reahzed decades ago that 
m simple 1-dimensional collapse models the outward ra- 
diation force on the accreting material should be signifi- 
cantly stronger than the inward pull of gravity ( Larson fc 



Starrfield|197l||Kahn|1974l|Wolfire fc (JassinelH|li)g7| m 
particular when taking dust opacities into account, ti mce 
stars as massive as 100—150 Mf?^ have been obse rved ( Bo- 
nanos et al.|2004[|Figer|2005[|Rauw et al.|2005| a sirhpe 
spherically symmetric approach to high-mass star forma- 
tion is doomed to fail. 

Consequently, two different models for massive star 
formation have been proposed. The first one takes ad- 
vantage of the fact that high-mass stars always form as 
members of stellar clusters. If the central density in the 
cluster is high enough, there is a chance that low-mass 
protostars co llide and so successiv ely build up more mas- 
sive objects (Bonnell et al. 1998| . As the radii of pro- 
tostars usually are considerably larger than the radii of 
main sequence stars in the same mass range this could 
be a viable option. However, the stellar densities re- 
quired to produce very massive stars are still extremely 
high and seem inconsistent with the observed values of 



Galactic star clusters (e.g. Portegies Zwart et al. 2010 



and references therein). An alternative approach is to 
argue that high-mass stars form just like low-mass stars 
by accretion of ambient gas that goes through a rota- 
tionally supp orted disk caused by angular momentum 
conservation (McKee & Tan 2003). Indeed such disk 
structures are observed around a number of high-mass 



Davies et a 



protostars (Chini e t al. 
|i>010| 



2004 2006 [Jiang et al. 2008 



Dreaks any spheri 



Their presence 
cal symmetry that might have been present in the initial 
cloud and thus solves the opacity problem. Radiation 
tends to break out along the polar axis, while matter 
is transported inwards through parts of the equatorial 
plane s hielded by the disk. Hydrodynamic s imulations 
in two (lYorke fc Sonnhalterl2002J and three ( ]KrumhoIz 
et al. 2009, Kuiper et al. 2010 in prep.) dimensions 



using a fiux-limited dittusion approach to the transport 
of non-ionizing r adiati on strongly support this picture. 
Krumholz et al. (20091 find non-axisymmetric Rayleigh- 
I'aylor instabil i ties us ing gray radiation transfer, whereas 



Kuiper et al. (2010), using frequency-dependent radi- 
ation transfer, do not find such instabilities, but nev- 
ertheless find strong disk accretion. Strong accretion 
continues when ionizing radiation is in cluded using ray 



tracing methods in three dim ensions (Dale et al. 2005 



2007 Peters et al. 2010a|b). In a clustered environ 



ment or m individual collapsing cores where the disk 
becomes gravitationally unstable, the three-dimensional 
models show that material flows along dense, opaque fil- 
aments whereas the radiation escapes through optically 
thin channels of low-density material. It has been demon- 
strated that even ionized material can be accreted, if the 
accretion fiow is strong enough. H ii regions are gravita- 
tionally trapped at that stage, but soon begin to rapidly 
fluctuate between trapped and extended states, in agree- 



ment with ob servations (Peters et al. 2010a Galvan 



Madrid et al. 2010) . Over time, the same ultracompact 



H II region can expand anisotropically, contract again 
and take on any of the obs e rved morphological classes 



(Wood fc Churchwell 1989 Kurtz et al. 1994 



et al.|2010a|b[ ). In their extended phases, expanding H li 



Peters 



regions drive bipolar neu tral outflows characteristic of 
high-mass star formation ( Peters et al. [[2010a) . 

Another key fact that any theory of massive star forma- 
tion must account for is the apparent presence of an up- 
per mass limit. No s tar more mass ive than 100— 150 M© 
has been observed (Massey 2003). This holds for the 
Galactic field, however, it is also true for young star clus 
ters that are massive enough so that purel y random sam 
pling of the ini tial mass function (IMF) ([Kroupa 



2002 



Chabrier 2003) without upp er mass limit should have 
yielded stars above 150 M© fWeidner fc Kroupa' 12004 



Figer 2005; Ocy fc Clarke 2005; Weidner ct al. 201U1 see 
however, [Selman fc Melnick[[200 8). This immediately 
raises the question ot what is the physical origin of this 
apparent mass limit. It has been speculated before that 
radiative stella r feedback might be resp onsible for this 
limit (see, e.g., Zinnecker fc Yorke[[2007 ) or alternatively 
that the internal stability limit o f stars with non-zero 
metal licity lies in this ma ss regime ( Appenzeller|1970b|a 



1987; jBaraffe et alTj"2001 1. However, fragmentation coulc 
also limit protostellar mass growth. Indeed, this is what 
we see in the simulations discussed here. The likelihood 
of fragmentation to occur and the number of fragments 
to form depends sensitively on the physical conditions 
in the star-forming clou d and its initial a.nd en vironmen- 
tal parame ters (see, e.g., Krumholz et al. 2010 JGirichidis] 



et al. 2010). Understanding the build-up of massive stars, 
therefore, requires detailed knowledge about the physi- 
cal processes that initiate and regulate the formation and 
dynamic al evolution of the m olecular clouds these stars 
form in ( Vazquez-Semadeni et al. 2009) . 

We argue that ionizing radiation, just like its non- 
ionizing, lower-energy counterpart, cannot shut off the 
accretion flow onto massive stars. Instead it is the dy- 
namical processes in the gravitationally unstable accre- 
tion flow that inevitably occurs during the collapse of 
high-mass cloud cores that control the mass growth of 
individual protostars. Accretion onto the central star is 
shut off by the fragmentation of the disk and the for- 
mation of lower-mass companions that intercept inward 
moving material. We call this process fragmentation- 
induced starvation and argue that it occurs unavoidably 
in regions of high-mass star formation where the mass 
flow onto the disk exceeds the inward transport of mat- 
ter due to viscosity only and thus renders the disk unsta- 
ble to fragmentation. We speculate that fragmentation- 
induced starvation is important not only for present-day 
star formation but also in the primordial universe during 
the formation of metal-free Population HI stars. Conse- 
quently, we expect these stars to be in binary or small 
number multiple systems and to be of lower mass than 



usually inferred (Abel et al. 2002[ Bromm et al 



Indeed, current numerical simul ations provide tne first 



2009) 



hints that this might be the c ase ( Clark et al.|2008" Turk 
eraL[[2009l [Stacy et al.[[2010l ) 



In the current s tudy, we analyze the simulations by[Pe 



ters et al. (2010a) with special focus on the mass growtl 
history of the individual stars forming and the physical 
processes that influence their accretion rate. We briefly 
review the numerical method we use and the assump- 



tions and approximations behind it in Section [2] Then 
we describe our findings in Section [3] and discuss them in 
the context of present-day and primordial star formation 
in Section m We summarize and conclude in Section [Sl 



2. METHOD AND ASSUMPTIONS 



radiation- 
formation 



We present three-dimensional, 
hydrodynamical simulations of massive star 
that include heating by ionizing and non- ionizing 
radiation u sing the adaptive- mesh code FLASH ( Fryxell 
et al.||200"0 ). We use our improved version of the hybrid- 
chara!cteristics raytr acing method ( Rijkhorst et al.J [2006 



Peters et al. 2010a) to propaga te the radiation on th e 
grid and couple sink particles |Federrath et al. 2010), 
which we use as models of protostars, to the radiation 
module via a prestellar model ( ^Peters et al.^^2010a ). 

The simulations start with a luOO Mq molecular cloud. 
The cloud has a constant density core of p = 1.27x lO^^'' 
^ within a radius of r = 0.5 pc and then falls off 
1.6 pc. The initial temperature of the 



gem 

as r~^/^ until r 

cloud is T = 30 K. The whole cloud is set up in solid body 

rotation with an angular velocity w = 1.5 x 10"^^ s^^ 

corresponding to a ratio of rotational to gravitational 

energy /3 — 0.05 and a mean specific angular momentum 

of j = 1.27 X lO^^cm^s-i. 

We follow the gravitational collapse of the molecular 
cloud with the adaptive mesh until we reach a cell size 
of 98 AU. We create sink particles at a cut-off density of 
Pcrit = 7 X lO^^^gcm""^. All gas within the accretion 
radius of Tgink — 590 AU above pcrit is accreted to the 
sink particle if it is gravitationally bound to it. The Jeans 



mass on the highest refinement level is Mw 



0.13 M, 



0- 

The adaptive mesh technique allows us to resolve the 
gravitational collapse of the gas from the parsec scale 
down to a few hundred AU. At higher spatial resolution 
of only several ten AU, the gas becomes optically thick 
to non-ionizing radiation, and scattering effects must be 
taken into account. Since our cut-off density is more 
than two orders of magnitude smaller than the onset of 



Larson 



the o ptically thick regime at ^ 10~ "^ gem"'' (e.g 

1969 1 and we are focussing in our analysis on large-scale 
"effects on the stellar cluster scale, we expect the raytrac- 
ing approximation to be valid. Feedback by radiation 
pressure, which plays a role on the very small scales and 
is neglected in our model, is dynamica lly unimportant on 
these large scales (Krumholz & Matzncr 2009). 

We discuss three simulations (see I'ablefTl). In the first 
simulation (run A), we only allow for the formation of 
a single sink particle and suppress the formation of sec- 
ondary sink particles artificially by introducing the dy- 
namical temperature fioor 



T'min = —j—p{nAxf 
tt/cb 



(1) 



with Newton's constant G, mean molecular weight /i, 
Boltzmann's constant fee, local gas density p, and cell 
size Ax. The temperature fioor ensures the sufficient res- 
olution of the Jeans leng th with n > 4 cells, w hich avoids 
artificial fragmentation ( Truelove et al.|[l997 |. The sup- 
pression of secondary sink tormation guarantees that the 
accretion flow around the massive star is not weakened 
by the fragmentation of the disk, which would otherwise 
inevitably lead to the formation of companion stars that 



limit thegrowth of the massive stars in the cluster (see 
Section [sT. In the second simulation (run B), this dy- 
namical temperature floor is not applied, and a whole 
stellar cluster forms during the course of the simulation. 
The third simulation (run D) is a control run that al- 
lows us to study the influence of radiation feedback on 
the stellar cluster evolution. As in run B, a small stel- 
lar cluster forms, but the stars emit neither ionizing nor 
non-ionizing radiation. 

To demonstrate convergence of our study, we have also 
run an additional simulation at twice the numerical res- 
olution. This simulation, run D-I-, which we ran for 
0.66 Myr, is identical to run D, except that the cell size at 
highest grid refinement is 49 AU, the sink particle radius 
is Tsink = 312 AU and the sink particle cut-off density is 
Pcrit — 2.5 X 10~^^gcm^'^. Since we ran it for only a 
fraction of the time covered by the other simulations, we 
will not discuss it at length, but the model ran for long 
enough to show that the relation between total stellar 
cluster mass and the mass of the largest star is the same 
is in the lower resolution simulations (see Section 4.2.2 
for further discussion). 

The numerical method along wit h its inherent physica l 
limitations is discussed in detail in Peters et al. (2010a). 



3. RESULTS 

3.1. Accretion History 

In this section we compare the protostellar mass 
growth rates from our three runs, with a single sink 
particle (run A), multiple sinks and radiative heating 
(run B), and multiple sinks with no radiative heatin g 
(run D). As already discussed by Peters et al. (2010a), 
when only the central sink particle is allowed to form 
(run A), nothing stops the accretion flow to the cen- 
ter. Figure [l] shows that the central protostar grows at 
a rate M « 5.9 x 10~^ Mq yr"-'^ until we stop the cal- 
culation when the star has reached 72 Mq . The growing 
star ionizes the surrounding gas, raising it to high pres- 
sure. However this hot bubble soon breaks out above 
and below the disk plane, without affecting the gas flow 
in the dislP] midplane much. In particular, it cannot 
halt the accretion onto the central star. Similar findings 
have also been reported from simulations focussing on 
the effects of non-ionizing radiation acting on somewhat 
smaller scales ("York e fc Sonnhalter 2002 Krumholz et al. 
2007b, 2009; Sigalotti et al.||2009l). Radiation pressure is 



not able to stop accretion onto massive stars and is dy- 
namically unimportant, except maybe in the centers of 
dense star cluste rs near the Galactic center (Krumholz 
fc Matzner||2009| . 



I'he situation is different when the disk can fragment 
and form multiple sink particles. Initially the mass 
growth of the central protostar in runs B and D is com- 
parable to the one in run A. As soon as further proto- 
stars form in the gravitationally unstable disk, they begin 
to compete with the central object for accretion of disk 
material. Howeve r, unlike in the cla ssical competitive 



accretion picture (Bonnell et al. 2001a, 2004), it is not 



^ We will hereafter refer to the flattened, dense, accretion flow 
that forms in the midplane of our rotating core as a disk. However, 
it is not necessarily a true Keplerian, viscous, accretion disk, which 
probably only forms within the central few hundred astronomical 
units, unresolved by our models. 



TABLE 1 
Overview of collapse simulations. 



Name Resolution Radiative Feedback Multiple Sinks Mginj^g (Mq) Wsini^s Mmax (Mq) 



Run A 
Run B 
Run D 



98 AU 
98 AU 
98 AU 



yes 
yes 
no 



no 
yes 
yes 



72.13 
125.56 
151.43 



1 

25 
37 



72.13 
23.39 

14.64 



Note. — We have run a simulation with a single sink particle (run A) to probe the upper mass limit with radiative feedback as well 
as simulations with multiple sink particles to study stellar cluster formation with (run B) and without (run D) radiative feedback. The 
simulations differ largely in the total mass in sink particles (Mginitg), the total number of sink particles (JVa j pyg) an d the mass of the 
most massive star (Mmax) at the end of the simulation. We follow the notation introduced in [Peters et al.| | 2010a^ , where additional 
low-resolution runs Ca und Cb were discussed. 
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Fig. 1. — top Total accretion history of all sink particles com- 
bined forming in runs A, B, and D. While the heating by non- 
ionizing radiation does not affect the total star formation rate, the 
ionizing radiation appreciably reduces the total rate at which gas 
converts into stars once the most massive object has stopped ac- 
creting and its H II region can freely expand. This effect is not 
compensated by triggered star formation at the ionization front, 
which is never observed during the simulation. The slope of the 
total accretion history in run B goes down because the massive 
stars (dashed line) accrete at a decreased rate, while the low-mass 
stars (dotted line) keep accreting at the same rate, bottom Instan- 
taneous accretion rate as function of time of the first sink particle 
to form in the three runs. These are generally the most massive 
sink particles during most of the simulation. While the accretion 
onto the star in run A never stops, the massive stars in run B 
and D are finally starved of material. Since the radiative heating 
suppresses fragmentation in run B, the final mass of the star is 
almost twice as high as in run D. 



the most massive object that dominates and grows dis- 
proportionately fast. Figure [T] shows that, although the 
accretion rates of the most massive stars {M > 10 M©) 
steadily decrease, the low- mass stars (M < 10 M©), keep 
accreting at the same rate. Although the detailed mass 
distributio n of the low-mass stars depends on numerical 
resolution ( Federrath et al.|[2010 ) , the net effect of their 
accretion should not, so long as they can form at all. 
The successive formation of low-mass objects in the disk 
at increasing radii limits subsequent growth of the more 
massive objects in the inner disk. Material that moves 
inwards through the disk accretes preferentially onto the 

sinks at larger radii. 

The key to the process, as already discussed by |Bate| 
(20001, is that fragments that form later at larger radii 
will preferentially accrete larger- angular momentum ma- 
terial. This will happen on a much shorter time than the 
timescales required to redistribute this angular momen- 
tum to other parts of the disk at larger radii by viscous 
or gravitational torques. Similar effects are discussed in 



the disk fragmentation studies by Kratter et al. (2010). 
This behavior is found in models of low-mass proto- 
binary disks, where again the secondary accretes at a 
higher rate than the primary. Its orbit around the com- 
mon center of gravity scans larger radii and hence it en- 
counters material that moves inwards through the disk 
before the primary star. This dri ves the system towards 



equal mass es and circular orbits (Bate & Bonnell 1997 



Bate|20(J0). In our simulations, after a certam transition 



period hardly any gas makes it all the way to the center 
and the accretion rate of the first sink particle drops to 
almost zero. This is the essence of the fragmentation- 
induced starvation process. In run B, it prevents any 
star from reaching a mass larger than 25 M0 . The Jeans 
mass in run D is smaller than in run B because of the 
lack of accretion heating, and consequently the highest 
mass st ar in run D grows to less than ISM©. In com- 



parison, 



Krumholz et al. 



(2007b) described simulations 
starting with a ten times less massive core than ours, us- 
ing both an isothermal equation of state and including 
accretion heating. In both cases, objects with roughly 
half the mass of our most massive object formed. The 
isothermal case formed an only slightly less massive ob- 
ject than the one including accretion heating, just as in 
our models. The smaller final masses of their objects 
compared to ours seems mainly just to be a result of 
their smaller initial core masses. 

Inspection of Figure [l] reveals additional aspects of the 
process. We see that the total mass of the sink particle 
system increases at a faster rate in the multiple sink sim- 
ulations, run B and D, than in the single sink case, run 
A. This is understandable, because as more and more 
gas falls onto the disk it becomes more and more unsta- 



ble to fragmentation, so as time goes by additional sink 
particles form at larger and larger radii. Star formation 
occurs in a larger volume of the disk, and mass growth 
is not limited by the disk's ability to transport matter to 
its cent er b y gravitational or viscous torques (compare 
Section 3.2). As a result the overall star-formation rate 
is larger than in run A. 

Since the accretion heating raises the Jeans mass and 
length in run B, the total number of sink particles is 
higher in run D than in run B, and the stars in run D 
generall y reac h a lower mass than in run B (compare 



Section 4.2.2). These two effects cancel out to lead to 
the same overall star formation rate for some time. At 
one point in the evolution, however, also the total accre- 
tion rate of run B drops below that of run D. At time 
t sa 0.68 Myr the accretion flow around the most massive 
star has attenuated below the value required to trap the 
H II region. It is able to break out and affect a significant 
fraction of the disk area. A comparison with the mass 
growth of run D clearly shows that there is still enough 
gas available to continue constant cluster growth for an- 
other 50kyr or longer, but the gas can no longer collapse 
in run B. Instead, it is swept up in a shell surrounding 
the expanding H ii region. 

It is also notable that the expanding ionization front 
around the most massive stars does not trigger any sec- 
ondary star format ion, which suggests that triggered star 
formation (Elmegreen & Lada 1977) may not be as effi- 
cient as expected, at least on the scales considered here. 

Figure |2] shows the individual accretion histories of 
each of the sink particles in run B and run D. Radia- 
tive heating cannot prevent disk fragmentation but raises 
the local Jeans mass. Hence, as discussed above, much 
fewer stars form in run B than in run D, and the mass 
of the most massive stars in run B is higher. It is also 
evident from the figure that star formation is much more 
intermittent in the case with radiation feedback (run B). 
The reason for this behavior is that the star formation 
process is controlled by the local Jeans mass, which de- 
pends to a large degree on how the fil ame nts in the disk 
shield the rad iation (c ompare Section l3.2 [), as previously 
noted ' ^ ' ' "^ ~ 



by [ Dale et al.| ( [2005^ ,2007J and" Krumholz et al. 
. shielding can lower the Jeans rnass temporarily 



( |2009[ ). Shielding can lower the Jeans mass temporarily 
and thereby allow gravitational collapse that would not 
have occurred otherwise. 

The accretion histories also reveal large differences be- 
tween sink particles within the same simulation. For 
example, in run D a sink particle that forms around 
t w 0.688 Myr accretes very rapidly and becomes one of 
the most massive stars in the cluster while other stars 
only accrete sporadically, interrupted with long inter- 
vals of no accretion at all, and some others even accrete 
only shortly after their formation. There are several 
reasons for these differences. First, accretion can stop 
when stars are ejected from the cluster by iV-body inter- 
actions. This happens mostly to low-mass stars with 
masses M < 2Mq. Second, intermediate-mass stars 
can be starved of material by the formation of compan- 
ions. This is the essential element of the fragmentation- 
induced starvation scenario. A star that can grow to in- 
termediate masses is embedded in an accretion flow from 
a larger gas reservoir. If this reservoir is gravitationally 
unstable, it can collapse to form companions. This may 
explain why the intermediate-mass stars in run B seem 



to have shorter time intervals of no accretion. Since the 
Jeans mass is higher in this case, the formation of com- 
panions is suppressed, and this allows the star to accrete 
for longer times. Third, the stars are not fixed at a po- 
sition in the disk plane but move quickly within the disk 
by gravitational interaction with their surrounding dense 
gas and with nearby stars. Hence, they are not tied to 
the filament they form in but can move away from it. 
If the star moves into a low-density void in the disk, 
accretion will stop until the star is embedded in higher 
density gas again. The sink particle in run D mentioned 
above accretes so vigorously because it moves along with 
its parental filament over a long time and no companions 
form around it. Fourth, accretion can stop when the star 
resides within a low-density H ii region. This applies to 
the first two massive stars in run B. Their accretion flow 
had become so weak that the ionizing radiation was able 
to isolate the stars from the high-density gas in the disk 
and stopped accretion. 

3.2. Disk Growth and Disk Instability 

To quantify the accretion of mass onto the disk and the 
subsequent instability, we calculate the amount of mass 
in a control volume that encloses the disk at all times. 
This control volume is the same for all runs and has 
dimensions 0.24pcx0.24pcx0.015pc. We show in Fig- 
ure [3] the mass Mdisk of non-accreted gas contained in 
the control volume, the mass Mginks contained in all sink 
particles as well as the total mass Mtot = Afdisk + -^smks 
for runs A, B and D. The disk mass in run A is much 
larger than in the multiple sink runs B and D since the 
mass cannot be absorbed by secondary sink particles. 
Instead, the mass accumulates in the disk plane. Con- 
sequently, the star in run A is embedded in a strong 
accretion flow at all times, which facilitates accretion at 
an approximately constant rate despite radiation feed- 
back. The evolution in runs B and D differs from the 
mass growth in run A as soon as secondary sink particles 
form. The continuous formation of new sink particles 
in an expanding region around the central star in these 
runs keeps Mdisk almost constant. Figure [3] shows that at 
t « 0.6 4 Myr, outflow of m aterial from the central H 11 



region ( [Peters et al.| 2010a[ ) reduces Mdisk even further 
in run B, so that it falls slightly below the case without 
feedback (run D) at late times. However, accretion onto 
secondary sink particles had already cleared the central 
region, so that this is a relatively minor effect. The ac- 
cretion onto sink particles is not affected by the central 
ionization and stays at a constant rate until the H 11 
region can break free (compare Figure [I]) . 

We investigate the initial instability of the disk with 
an analysis of the Toomre Q-parameter 



Q 



2c^n 



ttEG 



with sound speed Cs, angular velocity f2, surface density 
S, and Newton's constant G. The disk is hnearly stable 
for Q > 1 and linearly unstable for Q < 1 ( Toomre|1964 
Goldreich & Lynden-Bell 19651. The initial phases o 



disk mstabihty m run B are displayed in Figure J4] The 
figure shows slices of density, temperature and Qin the 
disk plane for four different times. One can clearly see 
that the most unstable parts of the disk are the fllamen- 
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Fig. 2. — Individual accretion histories for run B and run D. The figure shows the stellar masses as function of time for all sink particles 
that form in run B {left) and run D (right). Because the Jeans mass is lower without radiative feedback (run D), many more sink particles 
form in run D than in run B, and the mass of the most massive stars is also lower. 

evidence for fragmentation and the presence o f substruc- 
ture on ~ 1000 AU scales (as also proposed by Krumholz 
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Fig. 3. — Disk growth and sink particle formation. The plot 
shows the mass Afjiak of non-accreted gas contained in a control 
volume around the disk, the mass M^\^^^ of all sink particles and 
the total mass Aftot = Afjigj; + M^\.^^^ for run A (green), run B 
(blue) and run D (red). The disk mass in runs B and D is kept 
almost constant by subsequent sink particle formation, while the 
disk in run A continuously grows. The deviation between the disk 
masses in runs B and D at late times is caused by ionization-driven 
outfiows, but these do not affect the total star formation rate. 

tary structures that form as the disk grows in mass. The 
heating by stellar radiation can suppress instability lo- 
cally, but shielding by the dense filaments prevents the 
whole disk from becoming stable and restricts the heat- 
ing to small regions near the center of the disk that are 
surrounded by filaments. This shielding makes it possible 
for star formation to progress radially outwards despite 
accretion heating by the stars. The disk remains suffi- 
ciently cool at the inner edge for gravitational instability 
to set in and star formation proceeds inside-out in the 
disk plane. Hence, the effect of the filamentary structures 
in the disk is twofold: they are so dense that they ren- 
der the disk unstable locally; and, because of their high 
density, they can effectively shield the radiation from the 
stellar cluster near the center of the disk, so that radia- 
tive heating does not stabilize the outer parts of the disk. 
Indeed, high-resolutio n observations of high-mass accre- 
tion disk candidates ( jBeuther et al. [2009 ) provide some 



elTaLl l pOOTal) ; 

Kmce the lilaments that shield the outer parts of the 
disk from radiation arc optically thick, with an optical 
depth for non-ionizing radiation of several tens, it is im- 
portant to estimate the degree to which our simulations, 
based on a raytracing technique to propagate the radia- 
tion on the grid, are affected by the lack of heating by 
diffuse radiation. To test this, we use the adaptive-mesh 
radiative transfer program RADMC-3E[5 RADMC-3D 
is based on the stan dard Monte Carlo method of Bjork^ 
man & Wood (2001 1 in combinati on with Lucy 's methoa 
of treating optically thin region s jLucy 1999 ). It is the 
succe ssor of the RADMC code (Duliemond & Uominikl 
20041 and has been used previously to generate maps 
of dus t emission from these simulations ( [Peters et aT] 
2010b| . We have calculated self-consistent dust temper- 



atures of the simulation snapshot shown in the last row 
of Figure [4[ Assuming that dust and gas temperatures 
in the disk plane are equal, we can compare the Monte 
Carlo dust temperature with the simulation gas temper- 
ature, as shown in Figure [5] The comparison demon- 
strates that direct heating aominates over diffuse radi- 
ation in regions more than '-^ 500 AU away from the 
stars. The diffusely heated regions lie completely within 
the region heated in any case by direct radiation, so our 
raytracing method accurately describes the shielding by 
the filaments of the cold disk region where secondary 
fragmentation proceeds. 

To illustrate the tendency of star formation to occur 
at increasingly larger disk radii, we show the disk ra- 
dius at which new sink particles form as a function of 
time for runs B and D in Figure [6[ Because the accre- 
tion heating sets in already with the very first stars that 
form in run B, sink formation is suppressed at small disk 
radii initially. The massive stars slowly spiral outwards 
with time, so that at i ~ 0.67 Myr their radiation can 
be shielded by filaments in the disk. Within these fila- 
ments, the gas then cools until local collapse sets in and 

^ http://www.mpia.de/homes/dullemon/radtrans/radmc- 

3d/index.html 
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Fig. 4. — Early phase of disk instability in run B. The panels show slices of density, temperature and the Toomre Q-parameter at four 
different times. Each frame shows the simulation time and the mass of the most massive star in the cluster. The black dots indicate the 
positions of sink particles. One can see how the stellar radiation initially heats up the disk locally, which enhances the stability. The dense 
filaments, however, shield the radiation, and the cold material within and behind the filaments becomes unstable again. 
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Fig. 5. — Comparison of heating by direct and diffuse radiation. 
The slices show the gas temperature from the simulation snap- 
shot in the last row of Figure [4| (top) and the Monte Carlo dust 
temperature generated by RAIJMC-3D for the same snapshot (bot- 
tom). The dust temperature from diffuse radiation lies below the 
gas temperature almost everywhere outside the region heated by 
direct radiation (contour). The black dots indicate the positions 
of sink particles. 

two sink particles form near the center of the disk. Once 
the filament has dissolved, the gas heats up again and no 
further sink particles form in the inner disk region. 

4. DISCUSSION 

4.1. Initial Conditions 

4.1.1. Density Profile 

As we argue above, the fragmentation behavior of the 
disk forming around the massive central star depends 
sensitively on the initial and boundary conditions, i. e. 
on the physical properties of the high-mass cloud core. 
One of the key parameters that determines whether frag- 
mentation becomes widespread during the collapse of a 
massive cloud core is its initial density profile. Numerical 
simulations indicate that density profiles with flat inner 
core are more susceptible to fragmentation, while cen- 
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Fig. 6. — Radius of sink formation as a function of time for run B 
(top) and run D (bottom). Because of the absence of radiation 
heating, both the Jeans length and the Jeans mass are smaller in 
run D than in run B, giving rise to the formation of more numerous, 
but generally lower-mass stars. In both simulations, sink formation 
gradually occurs at larger disk radii. The accretion heating by the 
first stars to form suppresses sink formation at small disk radii in 
run B until relatively late in the cluster evolution. 

trally concentrated cores (for example such as singular 
isothermal spheres with p oc r"^) usually form only one 
or at most a few objects ( Girichidis et al.||2010 ). 

The density structure ot prestellar cores is typically 
estimated through the analysis of dust emission or ab- 
sorption using near-IR extinction mapping of background 
starlight, mapping of millimeter/submillimeter dust con- 
tinuum emission, and mapping of dust a bsorption against 
the b right mid-IR background emission ( Bergin & Tafalla 



2007| ) . A main characteristic of the density profiles de- 
rived with the above techniques is that they require a 
central flattening. The density of low-mass cores is al- 
most constant within radii smaller than 2500 — 5000 AU 
with typical central densities of 10^ 



10^ cm"'^ (Motte 



et al.||i998[ [Ward-Thompson et al.||1999[ ). A popu 
proach is to describe these cores as truncated isotl 



ar ap- 
isothermal 
(Bonnor-Ebert) sphere (Ebert||1955{ lBong or|[T956 ) , that 



often provides a good fit to the data ( Bacmann et al. 
2001| |Alves et al.||200l"| [Kandori et al.pOOS; ). BoimoF 
Ebert spheres are equilibrium solutions of sell-gravitating 
gas bounded by external pressure. However, such den- 
sity structure is not unique. Numerical calculations 



of the dynamical evolution of supersonically turbulent 
clouds show that transient cores forming at the stagna- 
tion points of c onverge nt flows exhibit similar morphol- 



OEL 



ffiallesteros-Parcdes et al.l2003 Klessen et al. 2005 



Banerjee et al. 20091). 



i'he situation is less clear when 

be- 



it comes to higli-mass cores (Beuther et al. 2007) 
cause most high-mass cores studied to date siiow at least 
some sign of star formation or turn out to consist of sev- 
eral sub-conde nsations when observed with high enough 
resolution (e.g. Beuther et al.||2005 Beuther & Henning 
[2009) . However, large- scale surveys, e.g. a s conducted m 
the Cygnus-X region (Motte et al. 2007), indicate that 



high-mass cores are in many aspects similar to scaled- 
up versions of low-mass cores. We followed these lines 
of reasoning and chose an initial density profile with flat 
inner core and r~^'^ density profile outside a radius of 
r = 0.5 pc. Although we begin with smooth rather than 
turbulent initial conditions, by the time ionizing radi- 
ation begins to be emitted, gravitational fragmentation 
has already produced substantial density perturbations. 

4.1.2. Rotation 

The second important parameter is the initial rotation 
of this core, which defines the total amount of angular 
momentum in the system and hence the radial extent of 
the protostellar accretion disk. To assess the relevance of 
fragmentation-induced starvation for high-mass star for- 
mation we therefore must compare our choice of /3 = 0.05 
for the ratio between rotational and potential energy in 
our 1000 AIq core both with observational data of molec- 
ular cloud cores and with numerical simulations of core 
formation. Our choice of values correspond to a specific 
angular momentum of j = 1.27 x lO^'^cm^s"^. 

The hig h-mass cloud cores observed by |Pirogov| 
|et al.| (J2003) with masses up to a few thousand solar 
masses show specific angular momenta j from a few 
xlO^^cm^s^^ up to lO^'^cm^s"^, corresponding to /3- 
values up to 0.07. We took this sample as motivation 
for our choice of parameters. Interesting in this con- 
text is also the large-scale rotational motion of molecu- 
lar gas around a clust er of hypercompact H li regions in 
G20.08-0.14 N seen by [Galvan-Madrid et al.|p009J . This 
flow has a radius of about 0.5 pc, roughly consistent with 
our sim ulation. Similar o bservations were rn ade by |Keto| 
( p90| for GlO.6-0.4 and [Welch et"all pMT] ) for WdMT 

It needs to be noted, however, that although the spe- 
cific angular momenta of observed cores exhibit consid- 
erable scatter there is a clear trend of decreasing j with 
decreasing core mass M. For example, the mean val- 
ues in the sample discussed by Goodman et al. ( 1993 1 
are j sa 6 x 10^^ cm^s~^ with masses up to several hun- 
dred solar masses, while Caselli et al. (2002) studied low- 
mass cores with typically only a tew solar masses and 
find j « 7 X lO^^cm^s""'^. In all cases the typical values 
for the ratio between rotational energy and potential en- 
ergy is /3 < 0.05. This trend continues down to stellar 
scales. A binary G star with a orbital period of 3 days 
has j « 10^® cm^ s~^, while the spin of a typical T Tauri 
star is a few xlO^'^cm^s"^. Our own Sun rotates only 
with j w lO^^cm^s"^. That means, during the process 
of star formation most of the initial angular momentum i s 
removed from the collapsed object (Bodenheimer 1995). 

If we look for guidance from numerical simulations of 
molecular cloud formation and fragmentation (^Mac Low[ 



& Klessen 2004[ Klessen et al. 2009) we see a very sim- 
iiar trend. 



Gammie et al. 



ind a typical mean 
specif ic angular mo mentum of j = 4 x lO^^cm^s"^, 
while Jappsen & Klessen ( 2004[ ) find values in the range 
10^"cm^s-^< 7<10" ' •^^^" ■ • ' 



^^ cm-^s ^ for low-mass cores in the 



siniulations by Klessen et al. ( 1998 ) and Klessen & Burk- 



ert 



(2001 2000p depending on the evolutionary s t ate of 



prot ostellar collapse. Se e also Tilley fc Pudritz[ ( [2007 1 
and [Offner et al. (2008). Altogether, we hnd that our 



initial conditions are consistent with simulations of core 
formation. 

4.2. Relation to Other Theoretical Models 
4.2.1. Monolithic Collapse Models 

The monolithic collapse model rests on the similarity 
between t he shape of th e observed core mass di s tribu- 



tion (M otte et al.| [T998[ [Johnstone et al.|[2000 



2006 



Lada "et al. 2008) a nd the stellar in itial mass function, 
liVll'' ( [Kroupa||2002j [Chabrier[[2003| . It assumes a one- 
to-one relilioiiDetweenTEe'HIstaButions with only a con- 
stant efficiency factor separating the two functions. Each 
molecular cloud core collapses to form a single star or 
at most a close binary system, with protostellar feed- 
back processes determini ng the efficiency of the process 
( Matzner fc McKee|2000 ) . This model reduces the prob- 
lem ot the origin of the IMF to the problem of determin- 
ing the mass spectrum of bound cores, although strictly 
speaking the idea that the IMF is set by the mass spec- 
trum of cores is independent of any particular model for 
the origin of that mass spectrum. Arguments to explain 
the core mass distribution generally rely on the statistical 



properties of turbulence (Klessen 2001; Padoan & Nord- 
lund[[2002[ [Hennebelle fc Chabrier||2008[ |200'9[ ), which 



generate structures with a pure powerlaw mass spectrum 
The thermal Jeans mass in the cloud then imposes the 
flattening and turn-down in the observed mass spectrum. 
However, there are a number of caveats. Many of the 
prestellar cores found in observational surveys appear to 
be stable entities and thus are unlikely to be in a state of 
active star formation. In addition, the simple interpreta- 
tion that one core forms on average one star, and that all 
cores contain the same numb er of thermal Jean s masses, 
leads to a timescale problem ( Clark et al.|[2007 ) that re- 
quires differences in the core rnass function and the IMF. 
Other concern about this model comes from hydrody- 
namic simulations, which seem to indicate that massive 
cores should fragment i nto many stars rather than col- 
lapsing rnonolithically (jDobbs et al.ll2005[ Clark &: Bon- 



et al. 



neli[[2006[ [Bonnell fc Ba te 2006; Fcdcrra tiTet al.[|2010[ ). 
'I'his objection is, however, weakened by the fact that 
magnetic fields are able to reduce the level of f ragmen- 
tation on sc ales of molecular clouds as a whole ('Heitsch 
as well as of collapsing cloud co res {tien- 



nebel 



|200T1) 



romang 



2008 



In addition, radiative feed 



Hennebelle fc Teyssi er 2008 T 
ack also is able to reduce the 



number of frag ments that form during the collapse of 
high-mass c ores (Krumholz et al. [[2007b I as well as low- 
mass cores ( Oftner et al. 2009| ) . But again, all current 
simulations indicate that neither accretion heating nor 
ionizing radiation can prevent the fragmentation of the 
massive dense disk that builds up during protostellar col- 
lapse, nor can t hey stop accretion onto the massive star 
([Yorke fc Sonnhalter 2002[ iKrumholz et al.[[2009| fPe'- 



10 



ters et aL]|2010a ). Instead, the heating merely increases 
the average mass of the fragments. This is also sup- 
ported by analytic estimates comparing typical accretion 
rates onto the disk with its ability to transport matter 
inwards thr ough viscou s and gravitat ional torques (Krat-[ 
ter fc Matzner 2006^ Kratter et al.p010| ). Our current 
results agree with these studies. 

4.2.2. Competitive Accretion 

A second model for the origin of the IMF, called com- 
petitive accretion, focuses on the interaction between 
protostars, and between a protostellar popula t ion and 
the gas cloud around it (Bonnc U et al.|[2001a|b Bonnell 
' fc Bate|[2002{ [Bate fc Bonnch 20n5| . in the competitive 
accretion picture the origin ot the peak in the IMF is 
similar to the monolithic collapse model, it is set by the 
Jeans mass in the prestellar gas cloud. However, rather 
than fragmentation in the gas phase producing a spec- 
trum of core masses, each of which collapses down to 
a single star or star system, in the competitive accre- 
tion model all gas fragments down to roughly the Jeans 
mass. Prompt fragmentation therefore creates a mass 
function that lacks the powerlaw tail at high masses that 
we observe in the stellar mass function. This part of 
the distribution forms via a second phase in which Jeans 
mass-protostars compete for gas in the center of a dense 
cluster. The cluster potential channels mass towards the 
center, so stars that remain in the center grow to large 
masses, while those that are ejected from the c luster cen- 
ter by JV-body interactions remain low mass (Klessen &: 
Burkert ||2000| [Bonneh et al.||2004[). In this model, the 
apparent similarity between the core and stellar mass 
functions is an illusion, because the observed cores do 
not correspond to gr avitationally bound structures that 
will c ollapse to stars ( Clark fc Bonnell|2006 Smith et al. 
^008| . 

The competitive accretion picture has been challenged, 
on the grounds that the kinematic structure observed 
in star-forming regions is inconsistent with the idea 
that protostars have time to interact with one another 
strongly before they completely accrete th eir parent cores 
( [Krumiiolz et al.||2005[ [Andre et aL]|2007[ ). 

'i'aken at face value, competitive accretion models show 
a correlation between the mass of the most massive star 
A^max and the total cluster mass Mginks during the whole 

2/3 

cluster evolution that is roughly M^ax oc MJ^^^ (Bon- 
nell et al. 20041. This correlation has been argued to 



represent a way to observationally co nfirm competitive 
accretion (Krumholz & Bonnell 200J7) and is in fact in 
good agreement with obse rvations ( [VVeidner fc Kroupa 



2006 Weidner et al. 



20101). However, we hnd that our 
simulations also reproduce the observed relation between 

Mmax and Mginks- 

Figure |7] shows Afginks as function of M^ax for run A, 
run B andrun D, the higher resolution convergence study 

2/3 

run D-l- and the relation M^ax = O.SOMgj^^j^g, which was 
found by Bonnell et al. (2004) as a fit to their simula- 
tion data. Over the whole cluster evolution, the curve for 
run D lies above this fit, while the curve for run B always 
lies below it. The fit agree s with our simulation data 
as well as it does to that of Bonnell et al. (2004). The 
higher resolution simulation run U-f was only carried out 



but its data also agrees very well with the fit, demon- 
strating convergence of the results from the lower reso- 
lution simulations. This indicates that the scaling is not 
unique to competitive accretion, but can also be found 
with the fragmentation-induced starvation scenario and 
hence cannot be used as an observational confirmation 
of competitive accretion models. 

The agreement of both models to this prediction is 
surprising because the accretion behavior of the most 
massive stars is totally different. Whereas the massive 
stars in competitive accretion simulations automatically 
accrete large fractions of the available gas because they 
reside in the center of the gravitational potential dur- 
ing the whole cluster evolution, the massive stars in 
fragmentation-induced starvation models have continu- 
ously decreasing accretion rates since they are starved of 
material by other cluster members and, in the final phase, 
by feedback from ionizing radiation. It seems that the 
observed relation between Mginks and Mmax is a very gen- 
eral result of protostellar interaction in a common cluster 
environment, and not an unambiguous sign of competi- 
tive accretion at work. 

We can gain further insight from Figure [7] It shows 
that the accretion heating suppresses low-mass star 
formation and that for all times the simulation with 
feedback contains a more massive star relative to the 
whole cluster mass than the simulation without feed- 
back. When the most massive star in run D reaches 
10 M0, much more gas is used up to create additional 
low-mass stars than is tunneled towards the most mas- 
sive one. Thus, the growth of the most massive star 
in run D is much more ineffective than in run B. Re- 



markably, all four simulations in Bonnell et al. (2004) 
show the turn-off towards accretion onto low-mass stars 
around 10 Mq, and only one simulation has formed a star 
more massive than 10 Mq. A very similar turn-off (indi- 
cated by an arrow in Figure u\ at the same mass scale 
is found in run D, but the radiative heating in run B 
shifts the turn-off towards higher masses, so that the 
scaling is followed closely over a larger mass range. In 
fact, since the turn-off occurs only at the end of the simu- 
lation and some of the massive stars are still accreting at 
this point, it is unknown up to what mass the scaling will 
continue. However, these results clearly indicate that ra- 
diative feedback is necessary to re produce the sca l ing fo r 
masses beyond 10 M© observed by Weidner et al. (2010). 



4.3. Relevance for Primordial Star Formation 

We note that the initial conditions adopted here may 
also be appropriate for the formation of metal-free Pop- 
ulation III stars. The very first generation of stars in the 
universe is thought to form at a redshift of z ~ 15 in rel- 
atively unperturbed and quiescent halos with masses of 
~ 10^ Mf7) in dark matter and about lO'^ Mf;^ in baryons 



( [Bromm fc Larson||2004[ [Bromm et al.||2009[ ). This gas 
has primordial coniposition, i.e. basically consists of hy- 
drog en and helium, and cools mostly via H2 line emission 



e.g. Anninos et al.|1997 Glover 2005 Glover fc Jappsen 



20071). I'hese halos ha ve a well-detined d ensity structure 
and rotational profile ( Abel et al. 2002 ) which leads to 
the build up of an exte nded accretion disk around the 
central object (see , e.g., Yoshida et al.||2008 for numeri 



cal simulations, or i'an fc McKee||20U4| tor an analytical 



for a fraction of the total time of the other simulations, model). Just as discussed above, if the mass flow onto 
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Fig. 7. — The total mass in sink particles Msinj^s as a function 
of the most massive star in the cluster Mmax- We plot the curves 
for run A, run B and run D, the convergence study r un D+ as well 
as the fit from the competitive accretion simulations [Bonneil et aT] 
([2004). All simulations follow the competitive accretion prediction 
with good accuracy. The maximum mass Mmax in run B is always 
larger than in run D for a fixed cluster mass A/ginks. The turn- 
off away from the scaling relation (indicated by arrows) is shifted 
towards higher masses by radiative feedback. The higher resolution 
study run D-|- was not followed long enough to show the turn-off. 

the disk exceeds its capacity for transporting material 
inward by viscous torques, the disk becomes suscepti- 
ble to gravitational instability and fragments. There are 
first indications that this indeed happened during Pop 
III star formation from high-resolution numerical simu 
lations that follow gravitational collapse bey ond the for 



mation of the central hydrostatic protosta r ( Clark et al. 
20081 iTurk et al.||2009l IStacy et al.||2010p. We propose 



that fragmentation-indiiced starvation may be the phys 
ical process that determines the masses of Pop III stars. 
T he situation is very different in atomic cooling ha- 



los ( |Wise fc Abel||2007| |Greif et al.||2009| . Once the 
virial mass of a halo exceeds a value of ~ 10* Mq , the in- 
falling gas gets partially ionized in the virialization shock. 
The enhanced abundance of free electrons triggers rapid 
II2 formation and some fraction of the gas cools down 
very rapidly. As a result there are streams of co ld gas 



that can tr avel all the way to the center (see also Dekel 
et al.||2009 ) and the overall velocity structure of the gas 
becoines highly chaotic and turbulent. We expect that 
these environmental conditions prevent the build-up of a 
large-scale accretion disk in the center of the halo and 
our model of fragmentation-induced starvation should 
therefore not be applicable. However, a more thorough 
analysis awaits the execution of detailed, high-resolution 
numerical simulations. 

5. SUMMARY AND CONCLUSION 

We have presented a detailed analysis of the 
fragmentation-ind uced starvat i on sce nario, which was 
first introduced by Peters et al. (2010a). We have studied 
the accretion history and disk instability in collapse sim- 
ulations of massive star formation. We have compared 
the results from a full simulation with radiation feedback 
and multiple sink particles with two control runs, one 
with multiple sink particles but without radiation feed- 
back and one with radiation feedback but only a single 
sink particle. The combined analysis of all three simula- 



tions allows us to establish fragmentation-induced star- 
vation as a workable model of massive star formation. We 
have also compared the new model to the monolithic col- 
lapse and competitive accretion models and speculated 
on the relevance of fragmentation-induced starvation for 
primordial star formation. 

The basic principle of fragmentation-induced starva- 
tion is star formation in a rotating gravitationally unsta- 
ble flow. Gravitational instability in such a flow leads 
to fragmentation of the accretion flow and the forma- 
tion of companions around the central massive star that 
starve the central star of infalling material. In our simu- 
lations, star formation occurs within a rotationally flat- 
tened, disk-like structure. As more and more material 
falls onto the disk, star formation proceeds radially out- 
wards, keeping the total disk mass roughly constant. The 
accretion history of individual stars is tightly connected 
with the position of these stars relative to the dense fil- 
aments and low-density voids in the disk. The fact that 
most if not all high-mass stars are observed in higher- 
order multiple stellar systems that themselves belong to 
more massive star clusters provides strong evidence for 
the widespread occurrence of fragmentation that forms 
the basis of our discussion. 

We find that accretion heating does not prevent frag- 
mentation of the disk, but leads to a higher local Jeans 
mass. As a result, fewer stars form with than without ra- 
diation feedback. The accretion heating shifts the masses 
of the most massive stars up, but leaves the average stel- 
lar mass almost unaffected. The accretion heating does 
not change the overall star formation rate of the whole 
stellar cluster. Feedback by ionizing radiation is unable 
to stop protostellar growth if the accretion flow is strong 
enough. However, if the accretion flow weakens due to 
starvation, an H 11 region can expand and terminate the 
accretion process. The growing H 11 regions reduce the 
total star formation rate and do not trigger star forma- 
tion at the ionization front. 

Ou r model is able to e xplain the observed m orpholo- 
gies (Peters et al. 2010b|) and time variability (Galvan- 
Madrid et al. 2010) of ultracompact H 11 regions. We 
hnd that we can consistently reproduce the observed re- 
lation between the total mass of the star cluster Mginks 
and the maximum stellar mass in the cluster Mmax, 

2/3 
Miuax oc MJ^^^. This relation seems to be the general 

outcome of protostellar interaction in a common cluster 
environment rather than being a signpost of competitive 
accretion only, as previously claimed. In fact, the dynam- 
ical processes discussed here exhibit exactly the opposite 
behavior of competitive accretion, rather than run-away 
accretion onto the most massive star together with the 
suppression of the growth of lower-mass objects, we see 
that angular momentum conservation and the presence 
of lower-mass objects limit the mass growth of massive 
stars. 

Our simulations provide evidence for the rejection of 
proposals that the observed maximum stellar mass of 
~ 100 Mq is set by radiative feedback. When disk frag- 
mentation is artificially suppressed (run A) the central 
protostar accretes material at very high rate unimpeded 
by the intense UV radiation it em its without any indi - 
cations of an upper limit (see also Peters et al. 2010a). 
When we permit disk fragmentation to occur, it is the 
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process of fragmentation-induced starvation that pre- 
vents the steUar mass to become larger than ~ 25 Mq 
with our choice of initial conditions. We expect more 
massive, more centrally-condensed, and/or slower rotat- 
ing cloud cores to lead to more m assive protostars. In- 
deed, extensive parameter studies ( Girichidis et al.|2010 | 
show that the initial density profile dominates the accre- 
tion behavior, explaini ng the formation of a 4 M© star 



from a 100 M© core in Krumholz et al. (2009). The aL 



ternative view is to attribute the apparent stellar mass 
limit to internal stability constraints. 
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